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Abstract. Consider n — 21 > A point particles with equal masses 
in space, subject to the following symmetry constraint: at each in- 
stant they form an orbit of the dihedral group Di , where Di is the 
group of order 21 generated by two rotations of angle tt around 
two secant lines in space meeting at an angle of tt/L By adding a 
homogeneous potential of degree —a for a e (0, 2) (which recovers 
the gravitational Newtonian potential) , one finds a special n-body 
problem with three degrees of freedom, which is a kind of gen- 
eralisation of Devaney isosceles problem, in which all orbits have 
zero angular momentum. In the paper we find all the central con- 
figurations and we compute the dimension of the stable/unstable 
manifolds. 
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1. Introduction 

The goal of this paper is to compute all the central configurations and 
the dimension of the stable/unstable manifolds for the dihedral sym- 
metric n-body problem in space under the action of a homogeneous 
potential of degree —a. For the Newtonian potential this problem is 
a kind of generalisation of Devaney planar isosceles three body prob- 
lem p!Ol [TT].The dihedral problem is a special case of the full n-body 
problem which reduces to a Hamiltonian system with three degrees of 
freedom. Briefiy, one takes n = 2/ > 4 equal masses whose initial po- 
sition and velocity are symmetric with respect to the dihedral group 
of rotations Di C SO (3). So the masses form a (possibly degenerate 
and non-regular) antiprism in space (and they are vertices of two sym- 
metric parallel /-gons). Because of the symmetry of the problem, the 
masses will remain in such a configuration for all time. Hence we have 
a system with only three degrees of freedom. For / = 2, the four bod- 
ies are at vertices of a tetrahedron, and the problem has been studied 
in a series of papers by Delgado and Vidal [2H [9]. The main tool is 
the use of McGehee coordinates introduced in [17] but for a general 
homogeneous potential of degree —a. We replace the singularity due 
to total collapse with an invariant immersed manifold in the full phase 
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space usually called total collision manifold which is the immersion 
of the parabolic manifold of the projected phase space. We explicitly 
compute all central configurations for this problem and show that just 
three types can arise: a planar regular 2/-gon, a regular /-gonal prism 
and a /-gonal anti-prism. 

The motivation in order to study this kind of problem is twofold. 
From one side this problem is difficult enough to put on evidence some 
chaotic behaviour of the full n-body problem and at the same time 
it is simple enough to carry out some explicit computations. From 
the other side the interest in this kind of problem is due to the fact 
that it includes a lot of other problems with two or three degrees of 
freedom studied in the past decades. The literature is quite broad and 
we limit ourself to quote only some closest results; among the others is 
the tetrahedral four body problem without and with rotation, studied 
respectively in [9J and ^24j , the rectangular four body problem studied 
by Simo and Lacomba in [21] . 

Acknowledgements. We are very grateful to the anonymous referees 
for their suggestions, comments and criticism which greatly improved 
the manuscript. 

2. McGehee coordinates, projections and regularisations 

Let y = M'' denote the Euclidean space of dimension d and n > 2 an 
integer. Let denote the origin G M*^. Let mi, . . . , m„ be n positive 
numbers (which can be thought as masses). The configuration space 
of n point particles with masses rrii respectively and center of mass 
in can be identified with the subspace of V"' consisting of all points 
q = (qi,...q„) G V"' such that J27=i^i^i ~ 0- ^et n denote the 
set {1, . . . , n} of the first n positive integers. For each pair of indexes 
i,jEn let Aj j denote the collision set of the z-th and j-th particles 
= {q ^ ^\Qi = Qj}- Let A = UjjAj j be the collision set. 

Let X C V"' be an open cone (MX = X) and let a > be a given 
positive real number. We consider the potential function (the opposite 
of the potential energy) defined by 




If M is the diagonal matrix, then Newton equations 




Mq = — 
oq 

1 form as 



(2.1) 
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where the Hamiltonian is H = H{q,p) = {-M ^p,p) — U{q). Then 

equations (2.1) can be written in polar coordinates by setting the mass 
norm in defined for every q G X as 

\\qf={Mq,q) 

and suitably rescaling the momentum as follows 

P = M\ 



s = — 

P 



p'p 



with a = 2(3. 



In these coordinates equations (2.1) can be read as 

(z,s)p 
s' = M'^z - {z,s)s 



' P' 



(2.2) 



I3{z,s)z + —{s), 



where the time has been rescaled by dt = p^^^dr (that is, — 

dr 

d 



P 



dr 



; now the energy can be written as 



:2.3) H=-p-^{M-'z,z)-p-^U{s) 



{^-{M-'z,z)-U{s)^ 



Let k := dn and let us consider the projection (g, p) ^ (s, z) from the 



full phase space X x M'^ to the reduced space S 
trivial M'^-bundle on the ellipsoid S^^^) 

X X M'^ ^ 5^-1 X 1 



k-l 



X 



(which is the 



^-^ X M^, that is 



In McGehee coordinates it is easy to see that the fiow on X x M'^ can 
be projected to S 

(2.4) 



{z, s)s 



= (3{z,s)z + ^{s). 

Also, being X a cone, it is a cone on its (/c — l)-dimensional intersection 
with the ellipsoid S^~^, which we will denote simply by S" = S^~^ flX. 
We define the parabolic manifold as the projection of all zero-energy 
orbits (or, equivalently, of the zero-energy submanifold of X x R*^) in 
S X R'^, that is 



P ■= {(s, z) G 5 X R^ : -{M-^z, z) = U{s)} C S' 



1 



X 



Its dimension is dim S + k — l = 2k — 2. This is also the projection of 
McGehee total collision manifold (see [171 Uni UHl EH]); the manifold of 
(s, z) here is not considered as embedded in the space of (p, s, z) with 
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p = 0. By the form of equation (2.2) , it is easy to prove the following 
proposition. 



2.1. Lemma. Solutions of (2.4 ) in S xM. are projections of solutions 



of (2.2). The parabolic manifold P is invariant for the flow of (2.4 ) 
and solutions in P can be lifted to X x M.'' by integrating the equation 
p'/p = {z,s). 

The parabolic manifold P is the boundary of the (2A; — l)-dimensional 

elliptic and hyperbolic manifolds, defined as 



Elliptic = {(s, z)eSxW : -{M'^z, z) < U{s)} C S' 



Hyperbolic = {(s, z) G x 



■.]^{hr'z,z)>u{s)}^s 



k-l 



X 



They are again invariant (even if the function ^^-^ ^) is not an 



invariant of the flow in 5* 



fe-i 



X 



and correspond to projection of el- 



liptic/hyperbolic orbits (that is, orbits with negative/positive energy). 
In fact, any fixed-energy (negative/positive) surface is homeomorphic 
to the elliptic/hyperbolic manifold. Given a solution of (2.4) in the 
elliptic or hyperbolic manifolds, for each energy value h the lifted so- 
lutions in X X M'^ can be found simply by applying (2.3) as 

tn ,^ _ {M-'z,z)-2Uis) 

(^•5) P - ^ • 

The parabolic manifold P is fiberwise homeomorphic to a trivial (k — l)- 
sphere bundle on S" C S^~^. 

The next change of coordinates, due to McGehee [17] (with a refer- 
ence to Sundman [22j), is needed for defining the Sundman-Lyapunov 
coordinate v and for the regularisation of the parabolic manifold P. 
Let to G 



X 



be defined by 



{z,s) 
M-^z 



{z,s)s. 



Then z = vMs + Mw and {w, Ms) = 0, and equations (2.4) can be 
replaced by 

|2 , 



(2.6) 




+ (3v' - aU{s) 



w\\s + {l3- l)vw + M-'WsU{s), 



where Vg denotes covariant derivative, i.e. the component of the gra- 
dient tangent to the inertia ellipsoid ||q|| = 1: 

dU 

VsU = —is) + aUis)Ms. 
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The parabolic manifold P is then defined by the equation 

v"^ + ||iuf = 2f/(s). 

The trivial bundle S X is simply decomposed as the sum of the 
normal bundle {s,v) of 5* in and the tangent bundle TS (with 



coordinates (s,iu)). By the first equation in (2.6) 



v' = \\wf + (3v^-aU{s) = (1-/5) ||icf + a (^{\\wf + v^) - U{s] 

can be deduced the well-known fact that for 0<Q;<2,fisa Lya- 
punov function on the fiow in the parabolic and hyperbolic manifolds, 
and therefore the fiow is dissipative (gradient-like). Moreover, the equi- 



librium points in (2.6) are the projections of the equilibrium points of 
(2.2) (and the projection is one-to-one in the parabolic manifold), which 
can be found as solutions of 

{f 2 = 2U{s) 
V.f/(s) = 
w = 0. 

Hence all equilibrium points belong to the parabolic manifold P. The 
constant solution in a central configuration s with f ^ = 2U{s) can be 
lifted to the full space as a homotetic parabolic orbit by integrating 
(back to the real time coordinate) 

, / , N l/(l+/3) 

p = ±p~^^m^ =^ p{t) = (±(1 + /?)v/2f7(i)tj 

assuming the total collision occurs at t = (the -|- sign yields an 
ejection solution, the — sign yields a collision solution). More generally, 
homotetic solutions (i.e. s' = 0, with s(t) = s) can be found in the 
hyperbolic and elliptic manifolds by setting in equations (2.6) s' = w = 
0, and therefore by integrating the single equation 

v' = (3v^ - aU{s) 

and then lifting the solution found to the full space using the energy 



relation (2.5) The graphs of homotetic solutions are straight lines 



contained in the normal bundle of S" in S" x 



3. The dihedral 2n-BODY problem 

Let M'^ = C X M be endowed with coordinates {z,y), z E C, y E M. 
For Z > 1, let denote the primitive root of unity Q = e^'^*/'; the 
dihedral group Di C 50(3) is the group of order 21 generated by the 
rotations 

Cr- iz,y) iCiz,y) and k: {z,y) t-> iz,-y), 

where 'z is the complex conjugate of z. The non-trivial elements of 
Di = {Ci,k) are the / — 1 rotations around the /-gonal axis C/, j = 



ON THE DIHEDRAL n-BODY PROBLEM 



6 





(a) 1 = 2 



(h) 1 = 3 



Figure 1. Dihedral groups Di, with the upper half of 
the fundamental domains in white. 



1, . . . , / — 1 and the / rotations of angle vr around the / digonal axes 
orthogonal to the /-gonal axis (see figure [T]) C/'^; J = 1; • • • 5 ^- In fig- 
ures 1(a) and 1(b) one can find the upper-halves of the fundamental 
domains for the action of Di restricted on the unit sphere. In fact, in 
figure 1(a) corresponding to the dihedral four body problem, the fun- 



damental domain is represented by an octant of the shape sphere while 



figure 1(b) represent the fundamental domain on the shape sphere for 
the dihedral six body problem. 

Consider the permutation representation of Di given by left multi- 
plication (that is, the Cayley immersion a: Di ^ Tj21 of Di into the 
symmetric group on the 21 elements of Di, defined by a{g){x) = gx 
for each g,x E Di, see [I3] for more details). The action of Di on M.^ 
induces an orthogonal action on the configuration space R^' of n = 21 
point particles qi G in the three-dimensional space. The Newto- 
nian potential for the ra-body problem, homogeneous with degree —a 
induces by restriction on the fixed subspace (R^') ' = R^ a homoge- 
neous potential defined for each g G R'^ by 



(3.1) 



3eA\{i} 



gqi 



provided we assume (without loss of generality) all masses mf = l/l. 
Now, the potential U in (3.1) can be re-written in terms of coordinates 
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q = {z,y) e C xR as 

i-i I 

3=1 3=1 

By definition, for each g E Di, U{gq) = U{q). Further symmetries of 
U are: 

(i) the refiection on the plane y = (given by h: {z, y) ^ {z, —y)), 
(a) the / refiections on the planes containing the /-gonal axis and 

one of the digonal axes, 
(in) and the / refiections on the planes containing the /-gonal axis 
and the points (C/e'^'/', 0), j = 1, . . . , / . 

It is not difficult to prove that these are (up to conjugacy and multi- 
plication with elements in Di) all the elements of the normaliser of Di 
in 0(3). Thus we can study U only in the left-upper area of the Di- 
fundamental domain on 5^ C M^, as we have seen in figures 1(a) and 



1(b) Now, in order to simplify the expression of the potential we intro- 



duce the variables r and ^ as follows. lfy>0 and 2; 7^ 0, let r = r{z, y) 



be defined as r = 1 + 2?/^/|zP — 2y /\z\\/y'^ /\z\'^ + 1 and ^ = -. Hence 

z 

r G (0, 1], with r = 1 if and only if y = 0, 1 + = r(2 + %^/|z^|) and 
therefore 



CM' + %^ = kr ( ii - + = ^ |i - c/r^r 

1 ^ \^\Z If' ' 



(3.2) U 



In these coordinates the potential function U (q) can be written as 

Eii-c/i-"+r-"/^Eii-c/<r 

-J=l J=l 

We can now state the integral representation of the potential (3.2) 
proven in the Appendix El (see also [2J and remark IA.5I below) . 

3.1. Proposition. For /3 g (0,1), r g (0,1] and ^ e c C the 
potential U can be written as 



U 



, _, sin(/3vr) f' jl-tyf^t^-^ l-{trr 
' vr io {\-tT-f |l-(tr)'ef 

h /"^l-a 



•where ci is the constant ci = 1 1 ~ C/ 



Proof. For the proof of this result, see Appendix 1X1 q.e.d. 

3.2. Remark. The above integral representation plays a fundamental 
role in order to find all the central configuration. In fact, otherwise the 
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expression of the potential given in formula (3.2) is quite difficult to 
deal with. 

3.1. Planar type central configurations. On the unit sphere 5* C 
(of equation \z\'^ + = 1), parametrised by {(p,6) e (— 7r/2,7r/2) x 
[0, 27r) with y = sin ip and z = cosyje*^, the (reduced to the 2-sphere) 
potential reads 
(3.3) 



u{e,ip)={2cosipy 



E(-fj .E(-^(^ 

j=l V / j=i V 



l+tan Lf 



and by Proposition 13.11 also as 



U{9,r) 



where I{r,6) 



2\ P 



4r 



, _3sin(/57r) 
ci + lr ^ — ^— ^ 



/(r, 



TT 



1 (1 - 



1 - itr) 



21 



lo {1-tr^f 1 + {tryi - 2{try cos{2W) 
with (just for ip G [0,n/2)) 

tany9 1 — sin (p 



dt. 



1 + 2 tan^ ip-2 



cos (fi 1 + sin (y9 



and hence 



1 — r 



sm Lf 



and cos^ ip 



4r 



1 + r 1 + 

In spherical coordinates, the symmetry reflections of U are (up to con- 
jugacy) 

(i) the reflection on the horizontal plane: h^: {9,ip) ^ {6, —ip), 
(a) the reflection on the plane containing the /-gonal axis and the 

digonal axis he: {9,ip) i— *• {—9,ip) 
(Hi) and the reflection on the plane containing the Z-gonal axis and 
the point {e'^'/\ 0), deflned as h'g : {9, Lp) ^ {tc/I- 9, ip). 

As direct consequence of the Palais' symmetric criticality principle, it 
follows that critical points of the restrictions of the reduced potential U 
to the 1-spheres of such flxed planes are critical points for the restriction 
of U to the sphere, and hence are central conflgurations for U . In fact, 
as already observed this 1-spheres are nothing but the spaces fixed by 
each of the refiections given in (z), (u) and (iii). In principle it can 
be exist other critical points for the restriction of the potential U to 
the sphere which do not lie in these fixed spaces. However if we are 
able to show that out of this 1-spheres the derivative of the potential 
is bounded away from zero, we have done. 

Now consider the derivative with respect to 9 of f/, which by Propo- 
sition 13.11 can be written as follows 

7r(2r)" ^ ' ' 



(3.4) 



^ = -4/2sin(2/^^)- 
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where I{r,6) is strictly positive and defined for {6,r) ^ (2/^7?//, 1), k 

dU 

integer. Hence for eacli r G (0, 1] the derivative is strictly negative 

TT IX IT JvTT 

for 9 G (0, — ) and strictly positive for 9 G (— , — ). It is zero for 9 = — 

(2k + l)7r 

and r G (0, 1) and 9 = — and r = 1. Thus, for = 0, we have 

proved the following proposition: 

3.3. Lemma (Planar 2/-gon). For any a G (0,2) central configurations 
which are h^p- symmetric are on the vertices (e'^^'^"'"-^)'^*/'^^'), 0) of the reg- 
ular 2l-gon. 

3.2. Prism type central configurations. Now we have to explore 
the cases 9 = kir/l and 9 = (2k + l)7r/(2/), which correspond respec- 
tively to prisms and antiprisms. The derivative of (3.3) with respect 
to if is 



dU ^ 213- ^^^'^ 



^ COS^{jTT/l — 9) 

{sm\jTT/l -9)+ tan2 ipY^\ 



^3 5^) dip (2cos(^)^ 

^ tan 09 r . / M 

dU 

The term in square brackets fei^p) has the same sign of — — , and since 

dip 

ci is a constant and each term of the sum is strictly monotone in (p, 
for each 9 the function feif) can vanish at most once in the interval 
(0, vr/2). Since the limit of the sum as — > vr/2 is zero and q is positive, 
there will be a unique zero in (0, vr/2) (for a fixed 9) for all the values 
9 such that limi^^o /6»(v^) < 0, i.e. 

hm V cosHjrr/l-9) ^ ^ g (sin^(,VO)-' • 

--'h (sin^(W^-^)+tanV)'^^ ^ 

Now, since lim<^^o foiv) = — oo, there exists a unique minimum ip for 
9 = kn/l, k = . . .21 — 1, corresponding to a prism. 

3.4. Lemma (Prisms). There are exactly 41 central configurations which 
are hg-symmetric (up to conjugacy), and they are precisely on the ver- 
tices of a prism: (cos (^'e'^'^*/', ± sin (p') . 

We observe that in the dihedral four body problem these kind of 
central configurations collapse to square type central configurations. 

3.3. Antiprism type central configurations. It is left to compute 
critical points for 9 = {2k + l)7r/(2/), that is, to find zeroes of fe{ip) for 

TT 

= — , or, equivalently, /i^-symmetric central configurations. 
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3.5. Lemma ( Antiprisms) . There are exactly 21 central configurations 
which are h'^- symmetric (up to conjugacy) and f 7^ 0. They are on the 
vertices of a prism: (cos (pe^^'''^^^'^'^/^^^\ ± sin 0) . 

We remark that in the four body problem the antiprism type central 
configurations reduce to tetrahedral type configurations. 

Proof. It suffices to show that 

' cos^Un/l-n/{2l)) ^ g ^,^^^j,/i^y^ . 



If ll/2\ denotes the greatest integer n < 1/2, that is 

L//2J 

then 



(/-l)/2 /odd 
1/2 I even, 



^ cos^(j7r//-7r/(2/)) cos^QV/ - 7r/(2/)) 

^ (sin^O'TT// - 7r/{2l))y'-' ^ (sin2(j7r// - vr/(2/)))^+^ ' 
On the other hand 

l-l ll/2\ 

J2 {sm\j7:/l))-' = 2 (sin2(W0)"^ + di, 
3=1 3=1 



where 



1 / even 
I odd 



Now then, since 

cos^ X 1 1 



[sm X) [sm x) (sm xj 

the conclusion would follow once we could prove that 

ll/2\ 

2j2Cj>di, 

3=1 

where 

1 1 

a 



' (sin2(j7r//-7r/(20))^+' {sin\j7r/l - 7r/{2l))f {siu^jt: /l)f' 

If Z = 2, it turns out that Ci = 2^ - 1 and hence 2 ^L'i^J Q = 2Ci > 
= ^2. If Z = 3, then Ci = 2"(3 - 3"^), which is greater than 2 
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Figure 2 . Central configurations in the upper-half fun- 
damental domain. 

for all a = 2/3, so that 2Ci > 1 = di. In general, since j < 1/2, 
sm{jn/l) > sin(j7r// — 7r/(2/)), and therefore 

1 2 



' {sm\jn/l - n/{2l))Y^' {sm\jn/l - n/{2l))Y 
_ 1 -2sin2(j7r//-7r/(2/)) 
" {sm\j7r/l-n/i2l))Y^'' 
The first term is estimated by 

l-2sin2(7r/(2/)) ^ vr^ 2-^ 

(sin2(vr/(2/)))^-^^ ^ 2P^ 7r"+2 



2// 



4/2 /2 

TT^ 4 



and all other terms Cj with j >2 are in any case greater than — 1 ; thus 
for / > 4 

U/2J 2 

E^^>C'i-(L2J-l)> ^1-2 + 1^1-2-^-^' 

and thus for all / > 4 we have 2 Xlj=i^ Cj > 2 > di, which concludes 
the proof. q.e.d. 



Since there are no other central configurations, by (3.4) 
summarise the results in the following proposition. 



we can 



3.6. Proposition. All central configurations in the dihedral 2n-body 
problem are symmetric for one of the three types of reflections h^p 
Lemma \3.3i kg Lemma \3.4\ or h'g Lemma \3.5[ They are represented 
in the (upper-half) fundamental domain on the sphere in figure\^ 
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In fact, in figure [2] it is drawn a geodesic triangle which represents 
the fundamental domain on the shape sphere for the dihedral n-body 
problem. In this figure are shown the three types of central configura- 
tions arising in the problem we are dealing with in the exact location 
together with. Moreover we observe that due to the symmetry con- 
straint only two types of collisions can occur. We denoted by the name 
/-adic collision and binary collision, meaning that in the first case two 
clusters of /-bodies simultaneously collide, while in the second case I 
clusters of 2 bodies simultaneously collide. This two types of collisions 
are all located on the same plane containing the planar central config- 
urations while the /-adic central configurations can be represented in 
the north and south pole of the shape sphere. 

Now consider equations (2.6) in coordinates {6, ip) on the sphere: 

ds ds 

we set Wi and W2 such that w = Wi— + , i-e. (since s = 

00 dip 

(cosv9e*^,siny9)). 



w = wi{i cos V9e*^, 0) + W2{— sin yje*^, cos p). 



Then 



and 



COS 



2 I 2 



VsU{s) 



1 dU ds dU ds 
cos^ p 86 do dp dp) 



Also, equations (2.6) become 



(3.6) 



Wn 



2 2 I 2,^2 

Wi COS p + W2 + pv 



Wi 
W2 



aUie,p) 



{(3 — l)vwi + 2tan(y9 W1W2 + 



dU 



COS' 

dU 



'p dO 



{p - l)vw2 - -Wi sin2v9 + 

2 dp 



The linearization at equilibrium points (central configurations) (2.7) 
represented by the 5x5 matrix L 



IS 



-2(3v 





























1 d^U 


d'^U 


cos^ p dO"^ 

d^u 


dOdp 
d^U 





dpdO 


dp'^ 





1 












1 



Thus the eigenvalues of the linearization can be computed in terms of 
the eigenvalues of the Hessian D'^U{s) of U{9,p): 
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3.7. Proposition. The eigenvalues of L, at a central configuration's 
(i.e. at the point {v,s, 0), where v = ±a/2?7 (s) ), are equal to the roots 
A of the equation 

+ (1 - (3)vX = 7 
for each 7 eigenvalue of the Hessian D'^U{s). 

By elementary calculations it follows from Proposition 13.71 that the 
equilibrium points (±A/2f/(s), s, 0) are hyperbolic when the Hessian 
D^U is non-singular at s, and that for each positive eigenvalue 7 > 
of D'^U there is a pair of real eigenvalues of L, Ai > 0, A2 < 0; for each 
negative eigenvalue 7 < of D'^U , there are two eigenvalues Ai, A2 of 
L with negative real part (Ai, A2 are real ii d = (1 — + 47 > 

and Ai = A2 if (i = 0). 



3.8. Proposition. All equilibrium points of (3.6) are hyperbolic. 



Proof. We just need to proof that the Hessian D^U is non-singular at 
s, if s is a central configuration. Since each central configuration s 
lies in the line fixed by a reflection (which is a symmetry of U), the 
matrix D^U is diagonal at s. So the result follows once we prove that 



(s)y^O^ -T—^(s). But by (3.4), since I(r,6) is strictly positive 



and regular in a neighbourhood of s, -;t^(s) 7^ 0. By (3.5), the same 
^ OU'^ ' ' 

holds for ——-(s). q.e.d. 

3.9. Proposition. The dimension of the stable (unstable) manifold of 
(l;, s, 0) with V = y^2U{s) > is 3 (2) if 's is a 2l-gon or a prism; it 
is 2 (3) ifs is an antiprism. The dimension of the stable (unstable) 
manifold of the point {—v, s, 0) with v = ^/2U{s) > is equal to the di- 
mension of the unstable (stable) manifold of {v,'s,0) . The intersection 
of the stable (unstable) manifold of {v,s, 0) with the parabolic manifold 
P has codimension 0(l)inPifv>0. It has codimension 1 (0) in P 
ifv<0. 

Proof. These facts follow directly from the stable/unstable manifold 
theorem and the above arguments on eigenvalues of L. The results are 
summarised in table [H q. e. d. 



Appendix A. An integral representation for U 

The aim of this section is to give a direct proof of the integral rep- 
resentation for the potential U used before in order to compute all the 
central configurations. 
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Table 





V 


CO 3 


c c 

CO S 

B B 
^ ^ 


2/-gon and prism 


> 
< 


3 2 
2 3 


3 1 
1 3 


anti-prism 


> 
< 


2 3 

3 2 


2 2 
2 2 



Dimensions of stable and unstable manifolds 



For / > 2, let P/ denote the /-adic Perron- Frobenius operator, defined 
on complex functions / : 5*^ C C ^ C by 



For each k E Z, 



3=0 



1/ 



^'^Z' if = mod / 
if A; ^ mod /. 



In terms of the /-adic Perron-Frobenius operator, the potential (3.2) 
can be written as 



1.2) 



U = \z\ 



where q is the constant q = X]j=i |1 ~ C/l " and P; (|1 — r^l ") (^') 

denotes the function P; (|1 — r^p") of argument ^ evaluated at In 

order to compute Pi (|1 — t^I "), we expand |1 — r^|~" in a double 
power series as follows. 

A.l. Lemma. For each r e (0, 1] and a = 2/3 > 

+cx> 
n=—oo 

with, for each n > 0, 



TT 
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Proof. 

|l-re|-"=(l-rO-^(l-rr^)-^ 

/ OO / ^\ 



h,k=0 



OO ^ 

E 



n=—oo 



\ 



-1)" E ("/) 

fc_ft=„ V / 



fe,/i>0 



bn 

Now, recall that for each (3 > and integer 

-(^\_( ^^n[n + P-i\_ r(iv + /3)r(i-/3) 
AT ^ ^ V ^ y r(iv + i)r(/3)r(i - /5) 

(-ir r(Ar + /j)r(l-/3) 

" r(/?)r(i-/?) ■ r(iv + i) 
= 5(1-/5,^ + /^) 

TT 

where B{x,y) denotes the beta function, defined as 
and wc have used the equalities 

r(/3)r(i - /3) 



sin(/37r) ' 



-P\ _ {-(3){-(3 - 1) ... - AT + 1) _ T{N + /5) 



AT 7 AT! ^ ' r{N+l)T{(3) 

and 

r(/3 + i) 



^A^y T{N + 1)T{13 - N + 1)' 
We can now use the integral representation of the binomial function 

(1.3) (]^^) = {-if'^^l\i-t)-^t^-H^dt, 
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h=0 ^ '"'Jo 



sin(/37r) 



h=o V "'0 V / 



,2/i 



sin(/?7r) 



•^^ h=0 O \ ' 



q.e.d. 



A. 2. Lemma. For eac/i /? e (0, 1), r G (0, 1] and integer I > 2 



+ 00 



Proof. It follows directly from equation (1.1) The convergence is easy 
to check. q.e.d. 

A. 3. Lemma. For each (3 G (0, 1) and r G (0, 1] and integer I > 2 



sin(/37r) {l-t)-H^-^ l-(tr) 



TT 



21 



{l-tr^f \l-{trYi\- 



dt. 



Proof. 



n=— oo n=0 n=l 



sin(/57r) 



r'" / (1 - t)-^t^-^t'"(l - tr")"^ dt ) C"" 



n=0 

oo 

n=l 

sin(/?7r) /•^(l-t)-^t'3~i 



TT 



TT 



(1 - tr2) 



/3 



?i=0 



The conclusion follows since 



71=1 



2« 



(it 



n=0 



n=l 



1 - (tr) 
|l-(trKr 



q.e.d. 



Thus we proved the following result. 
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A.4. Proposition. For (3 e (0,1), r e (0,1] and ^ e c C the 
potential U can he written as 



U 



ci + lr ^ / — — ,/ ^Y'.Lno dt 



sm{p7i) (1 - t)-H^'^ 1 - {trf^ 



IT 



(1-trY \l-{trye? 



where ci is the constant ci = Yl\=i 1 1 ^ C/ 



A. 5. Remark. An analogue of the integral representation of the poten- 
tial is well known, and can be traced back to F.-F. Tisserand's book [23] 
(chapter XVII) for the exponent a = 1; it had been used by M. Lin- 
dow [IBJ (section 3) in computing central configurations for the planar 
gravitational field generated by a regular n-gon. More recently D. Bang 
and B. Elmabsout extended and generalised Lindow's theorem, proving 
an equivalent of \3.1\ (Proposition 7 and 8 of [2\). The proof given here 
is direct, and allows explicit estimates that can be used to compute the 
Hessian for the potential restricted to the shape sphere. Furthermore, 
it involves an interesting connection with the l-adic Ruelle-Perron- 
Frohenius operator (see P. Gaspard's paper ^S\), which is worth a men- 
tion. 
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